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Abstract 

We show that the class of hyperelliptic solutions to the Ernst equation (the station- 
ary axisymmetric Einstein equations in vacuum) previously discovered by Korotkin 
and Neugebauer and Meinel can be derived via Riemann-Hilbert techniques. The 
present paper extends the discussion of the physical properties of these solutions that 
was begun in a Physical Review Letter, and supplies complete proofs. We identify a 
physically interesting subclass where the Ernst potential is everywhere regular except 
at a closed surface which might be identified with the surface of a body of revolution. 
The corresponding spacetimes are asymptotically flat and equatorially symmetric. This 
suggests that they could describe the exterior of an isolated body, for instance a rel- 
ativistic star or a galaxy. Within this class, one has the freedom to specify a real 
function and a set of complex parameters which can possibly be used to solve certain 
boundary value problems for the Ernst equation. The solutions can have ergoregions, 
a Minkowskian limit and an ultrarelativistic limit where the metric approaches the 
extreme Kerr solution. We give explicit formulae for the potential on the axis and in 
the equatorial plane where the expressions simplify. Special attention is paid to the 
simplest non-static solutions (which are of genus two) to which the rigidly rotating 
dust disk belongs. 



PACS numbers: 04.20. Jb, 02.10.Rn, 02.30.Jr 



1 Introduction 



It is generally believed that most of the stars and galaxies can be described in good ap- 
proximation as fluid bodies in thermodynamical equilibrium. In the framework of general 
relativity, this implies (see e.g. [p], 0) that the corresponding spacetimes are stationary and 
axisymmetric. Moreover it is usually assumed (though there is no proof known to us) that 
they are equatorially symmetric. This stresses the importance of the study of stationary ax- 
isymmetric spacetimes. A relativistic treatment is necessary for rapidly rotating and massive 
compact objects like pulsars, neutron stars and black-holes. 

Though the importance of global solutions describing stationary axisymmetric fluid bod- 
ies is generally accepted, the complicated structure of the Einstein equations with matter 
gives little hope that such solutions can be found in the near future. Only for special and 
somewhat unphysical equations of state ^ f|, [5], it was possible to give solutions in the mat- 
ter region which are discussed as candidates for an interior solution. In the exterior vacuum 
region, however, powerful solution generating techniques are at hand. Since the surface T z of 
a compact astrophysical object constitutes a natural boundary at which the metric functions 
are not continuously differentiable, one is looking for solutions to the vacuum equations that 
are analytic outside this contour and can be at least continuously extended to T z . This 
means that the typical problem one has to consider for the vacuum Einstein equations is a 
boundary value problem of Dirichlet, von Neumann or mixed type, see ||. The matter then 
enters only in form of boundary conditions for the vacuum equations. This is possible if an 
interior solution is known or if only two-dimensionally extended bodies like disks or shells 
are considered. In the latter case, the surfacelike distribution of the matter implies that the 
matter equations reduce to ordinary differential equations. Notice that disks are important 
models in astrophysics for certain types of galaxies. 

The reason why it is much more promising to treat only the vacuum case is the equivalence 
of the stationary axisymmetric Einstein equations to a single nonlinear differential equation 
for a complex potential, the so called Ernst equation @. The latter belongs to a family 
of completely integrable nonlinear equations that are studied as the integrability conditions 
for associated linear differential systems. The common feature of these linear systems is 
that they contain an additional variable, the so called spectral parameter, which reflects 
an underlying symmetry of the differential equations under investigation, in the case of the 
Ernst equation the Geroch group jS[ . Associated linear systems for the Ernst equations were 



given in || [10], |TT[. The existence of this parameter can be used to construct solutions by 
prescribing the singular structure of the matrix of the linear system with respect to the 
spectral parameter. 

One of the most successful solution techniques for nonlinear differential equations rests 
on methods of algebraic geometry and leads to the so called finite-gap solutions that can be 
expressed elegantly in terms of theta functions. Such methods were first used to construct 
periodic and quasiperiodic solutions to nonlinear evolution equations like the Korteweg- 
de Vries (KdV) and the Sine-Gordon (SG) equation. For a survey of this subject we refer 
the reader to |Tj| [13] . However, it was only recently that algebro-geometrical methods were 



applied to the Ernst equation, see ||14|| . The found solutions differ from similar solutions of 
other equations in several aspects, e.g. they are in general not periodic or quasi-periodic. 
The main difference is that this class is much richer than previously obtained ones. 
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The development of solution techniques yields a deeper insight into the structure of 
nonlinear differential equations. However, from a practical point of view, it would also be 
desirable to solve initial value problems or, for the Ernst equation, boundary value problems. 
One approach to solve boundary value problems of the above mentioned type with the 
help of the linear system is to translate the physical boundary conditions into a Riemann- 



Hilbert problem which is equivalent to a linear integral equation, see [12|. Neugebauer and 



Meinel |T5j succeeded in doing this in the case of the rigidly rotating dust disk. They were 
able to reduce the matrix problem on a sphere to a scalar Riemann-Hilbert problem on a 
hyperelliptic Riemann surface which can be solved explicitly via quadratures. By making 
use of the gauge transformations of the linear system we were able to show |T7j that this 
is possible in general if the boundary value problem leads to a Riemann-Hilbert problem 
with rational jump data. Up to now there is however no direct way to infer the jump data 
from the boundary value problem one wants to solve. The explicit form of the hyperelliptic 
solutions possibly offers a different approach to boundary value problems: one can try to 
identify the free parameters in the solutions, a real valued function and a set of complex 
parameters, the branch points of the hyperelliptic Riemann surface, from the problem one 
wants to solve. 

To this end we study a class of solutions - which is essentially equivalent to |T3| and 
jj| - that is constructed via a generalized Riemann-Hilbert problem on a hyperelliptic 



Riemann surface. We present a complete discussion of the singularity structure of these Ernst 
potentials. It is possible to identify a subclass of solutions that are everywhere regular except 
at some contour, which can possibly be related to the surface of an isolated body, where 
the Ernst potential is bounded. These solutions are asymptotically flat and equatorially 
symmetric, and thus show all the features one might expect from the exterior solution for 
an isolated relativistic ideal fluid. They can have a Minkowskian and an extreme relativistic 
limit in which the body is 'hidden' behind a horizon, and the exterior solution becomes the 
extreme Kerr solution. This provides the hope that further solutions to physically interesting 
boundary value problems to the Ernst equation, besides the rigidly rotating dust disk, can 



be identified within this class. First results on this subclass where published in [18 



The paper is organized as follows. In section we introduce the linear system associated 
to the Ernst equation and discuss how the matrix of the system has to be constructed in 
order to end up with new solutions to the Ernst equation. Using the results of [0, we 
show how Riemann surfaces arise naturally in the context of linear systems with a spectral 
parameter. In the case of the Ernst equation, these are hyperelliptic Riemann surfaces with a 
special structure of the branch points. We will restrict ourselves to regular compact Riemann 
surfaces and are eventually led to consider families of hyperelliptic Riemann surfaces of 
arbitrary genus, parametrized by the physical coordinates. 

In section 3 we recall some basic notions of the theory of Riemann surfaces, theta functions 
and the solution of Riemann-Hilbert problems on Riemann surfaces due to Zverovich, and 
present the class of solutions. It is shown that the solution of the axisymmetric Laplace 
equation which can be freely prescribed in []I6| is a period of the Abelian integrals which 
determine the singularity structure of the matrix of the linear system. The differential 
relations between these periods are a subset of the so called Picard-Fuchs equations which 
we write down for the Ernst equation. In section 4 we discuss the singularity structure of 
these solutions. It is shown that the solutions can have a regular axis and are in general 
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asymptotically flat. Using an identity for theta functions, we are able to give in section 
5 compact formulas for two metric functions and a simple condition for the occurrence of 
ergospheres. A subclass of solutions with equatorial symmetry is presented in section 6. The 
common physical features of this subclass like the extreme relativistic limit are discussed. 
In section 7, we use the equatorial symmetry to give simplified formulae for the potential 
in the equatorial plane and on the axis. Since the rigidly rotating dust disk belongs to the 
simplest non-static solutions which are of genus 2, we consider this case in detail in section 
8. In section 9, we summarize the results and add some concluding remarks. 



Linear System for the Ernst equation and Mon- 
odromy matrix 



It is well known (see ||20|| ) that the metric of stationary axisymmetric vacuum spacetimes 
can be written in the Weyl-Lewis-Papapetrou form 

ds 2 = -e 2U (dt + ad0) 2 + e~ 2U (e 2k (dp 2 + d( 2 ) + p 2 d0 2 ) (2.1) 

where p and ( are Weyl's canonical coordinates and dt and 8$ are the two commuting 
asymptotically timelike respectively spacelike Killing vectors. 

In this case the vacuum field equations are equivalent to the Ernst equation for the 
complex potential / where / = e 2U + ib, and where the real function b is related to the 
metric functions via 

b = -~e m a z . (2.2) 
P 

Here the complex variable z stands for z = p + i(. With these settings, the Ernst equation 
reads 

1 2 

fzz "I" 777 ; z^ifz "I" fz) = "7 ; jfzfz j (2-3) 
2(z + z) f + f 

where a bar denotes complex conjugation in C. With a solution /, the metric function U 
follows directly from the definition of the Ernst potential whereas a can be obtained from 
via quadratures. The metric function k can be calculated from the relation 

k z = 2p(U z ) 2 -^e iU (a z ) 2 . (2.4) 
2p 



The integrability condition of (2.2) and (2.4) is the Ernst equation. 

The remarkable feature of the Ernst equation is that it is completely integrable. This 
means that it can be considered as the integrability condition of an overdetermined linear 
differential system for a matrix valued function $ that contains an additional variable, the 
so called spectral parameter K. The occurrence of the linear system with a spectral param- 
eter is a consequence of the symmetry group of the Ernst equation, the Geroch group ||. 



Several forms of the linear system are known in the literature (P, ITU, IT]). They are related 



through gauge transformations (see WW). The choice of a specific form of the linear system 
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is equivalent to a gauge fixing. We will use the form of [11], 

^ z ' 2 ) = { ( o m ) + ( m o ) } $( ^' ^ z ' f) = w ® ' (2 - 5 a) 

= {(^^) + ^^y (j^)}$(^wM) = ^> ( 2 - 5b ) 

where 

^7T7> N = T+T (2 ' 6) 

Obviously M and iV depend only on the coordinates 2 and 2 and not on the spectral param- 
eter K that lives on the Riemann surface C(z, z) — C given by fJ^(K) = (K — iz)(K + iz). 
Notice that £ is a Riemann surface of genus zero with coordinate dependent branch points. 
This is a special feature of the family of chiral field equations to which the Ernst equa- 
tion belongs that has no counterpart among the completely integrable nonlinear evolution 
equations for which algebro-geometric solutions have been constructed first. 
On £ we have an involutive map a, defined by 



C 3 P = (K, ±^{K - iz)(K + iz)) -> a(P) = P CT = (K, T^J{K - iz){K + iz)) G C , 

(2.7) 

and an anti-holomorphic involution r, defined by 



C3 P=(K, ±^J(K - iz)(K + iz)) -> r(P) = P=(K, ±^{K - iz)(K + iz)) e C . (2.8) 

It is possible to use the existence of the above linear system for the construction of 
solutions to the Ernst equation. To this end one investigates the singularity structure of 
the matrices $ 2 $ _1 and with respect to the spectral parameter and infers a set of 

conditions for the matrix $ (at least twice differentiable with respect to z, z) that satisfies 
the linear system (|2.5 a|) and ( p.5 b|) . This is done (see e.g. [ |L"f 



in 



Theorem 2.1 Let $(P) (P £ £) be a 2 x 2 -matrix with the following properties: 

I. 3>(P) is holomorphic and invertible at the branch points Pq = —iz and Pq such that the 
logarithmic derivative $ 2 $ _1 diverges as (K + iz)* at P and $ f $ _1 as (K — iz)* at P . 

II. All singularities of $ on C (poles, essential singularities, zeros of the determinant of 
$, branch cuts and branch points) are regular which means that the logarithmic derivatives 
$ 2 $ -1 and are holomorphic in the neigbourhood of the singular points (this implies 
they have to be independent of z, z). In particular $(P) should have 

a) regular singularities at the points Ai £ C (i — 1, . . . ,n) which do not depend on z, z, 

b) regular essential singularities at the points Si (i — 1, . . . , m) which do not depend on z, z, 

c) boundary values at a set of (orientable, piecewise smooth) contours 1^ C C (i — 1, . . .,1) 
independent of z, z, which are related on both sides of the contours via 

<S>„(P) = $ + (P)g i (P)\ PeVi . (2.9) 

where Gi{P) are matrices independent of z, z with Holder-continuous components and non- 
vanishing determinant. 

III. $ satisfies the reduction condition 

$(P CT ) = <x 3 <l>(P)7 , (2.10) 
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where 03 is the third Pauli matrix, and where 7 is an invertible matrix independent of z and 
z. 

IV. The normalization and reality condition 

$(P = oc+)=({ ]). (2.11) 



Then the function f in ( \2.11j) is a solution to the Ernst equation. 



A proof of this Theorem may be obtained by comparing the above matrix $ with the linear 
system ( |2.5 a| ) and ( |2.5 b| ). 



Proof : Because of I, $ and $ 1 can be expanded in a series in t = yK + iz and t' = 
— iz in a neighbourhood of P = Pq and P = Pq 7^ Pq respectively at all points Pq, Pq 
which do not belong to the singularities given in II. This implies that $> z $>~ 1 = oto/t + ci\ + 
ait + • ■ ■. We recognize that, because of I and II, ^ z ^ x — ao/t is a holomorphic function. 
The normalization condition IV implies that this quantity is bounded at infinity. According 
to Liouvilles theorem, it is a constant. Since $, and $ 2 are single valued functions on C, 
they must be functions of K and [Iq. Therefore we have $ 2 $ _1 = (5qJ J^Epj*+/3i- The matrix 
/3q must be independent of K and /i since $ 2 $ _1 must have the same number of zeros and 
poles on C. The structure of the matrices (3$ and j3\ follows from III. From the normalization 
condition IV, it follows that $ 2 $ _1 has the structure of (|2.5 a ). The corresponding equation 



for $ 2 $ 1 can be obtained in the same way. □ 

For a given Ernst potential /, the matrix $ in the above theorem is not uniquely de- 
termined. This reflects the fact that the gauge is not uniquely fixed in the linear system 
( [2.5 a|) and ( [2.5 b|) . If we choose without loss of generality 7 = U\ (the first Pauli matrix), 
the remaining gauge freedom can be seen from 



Corollary 2.2 Let he a matrix subject to the conditions of Theorem \2. 1\ , and C(K) 

be a 2 x 1-matrix that only depends on K 6 C with the properties 

C(K) = + a 2 (.K>i, 

ai(oo) = 1, a 2 (oo) = 0. (2.12) 



Then the matrix $'(P) = Q(P)C(K) also satisfies the conditions of Theorem \2.1\ and 
$'(oo+) = $(oo+). 

It is this gauge freedom to which we refer when we speak of the gauge freedom of the linear 
system in the following. 

It is interesting to note that the metric function a can be obtained from a given matrix 



$ without solving the equation ( |2.2|) , see |TJ]. We get 



Proposition 2.3 Let 5 be a local parameter in the vicinity of 00 . Then 

(a - a )e 2U = i($ n - $ 12 ) >5 , (2.13) 

where a is a constant that is fixed by the condition that a = on the regular part of the 
axis and at spatial infinity, and where $^ denotes the linear term in the expansion of $ in 
S divided by 5. 



5 



The proof follows from the linear system ( [2.5 a| ) and ( 2.5 b ). 



Proof : It is straightforward to check the relation 

= . (2.14) 



With (|2.5 a[), we get 



from which, together with ( |2.2| ), (|2.13|) follows. □ 
Notice that a is not gauge independent (in the sense of the above corollary) whereas a 

is. 

Theorem 2A can be used to construct solutions to the Ernst equation by determining the 
structure and the singularities of $ in accordance with the conditions I-IV. For nonlinear 
evolution equations, large classes of solutions were constructed with the help of algebro- 
geometric methods, in particular Riemann surface techniques. A keypoint in this context is 
the occurrence of Riemann surfaces which are related to the linear system of the integrable 
equation under consideration. In this paper we want to show how solutions for the Ernst 
equation can be constructed by making use of the so called monodromy matrix of the Ernst 
system, which - following |19| - can be introduced as follows. 



For a given linear system ( |2.5 a| ) and ( |2.5 b| ), we define the monodromy matrix L as a 
solution to the system 

L Z =\W,L], L- Z = [V,L] . (2.16) 
For a known solution $ of ( 2.5 a|) and ( [2. 5 b|) , L can be directly constructed in the form 

L(K) = -fi^K^C®- 1 (2.17) 

where C is an arbitrary constant matrix with detC = — 1 and (x does not depend on the 
physical coordinates. Since $ is analytic in K, there is a solution to ( |2.16| ) with the same 
properties. 

It follows from ( |2.16|) that the coefficients of the characteristic polynomial Q(n, K) = 
det(L(K) — fil) are independent of the coordinates. Without loss of generality we may 
assume TtL(K) = 0. Then L has the structure 

T _(A{K) B{K)\ 

L -\C(K) -A(K) ) ■ (2 - 18) 

The equation Q(£l,K) = 0, i.e. 

fi 2 = A 2 + BC , (2.19) 

is then the equation of an algebraic curve which in general will have infinite genus. We will 
restrict the analysis in the following to the case of a regular curve with finite genus. 
In this case, the Riemann surface L is given by an equation of the form 
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A 2 = f[{K-E i )(K-F i ) (2.20) 
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El F?---E° F? 



iz 



-iz 



El Fl • • ■ Eg Fg 



Figure 1: The Hurwitz diagram of C 

where E± and Fi are obviously independent of the physical coordinates. This equation 
represents a two sheeted covering of the Riemann sphere and thus a four sheeted covering 
of the complex plane. A point P G C can be given by P = (K, fi (K), fi(K)). The Hurwitz 
diagram of L is shown in figure 1. 

There is an automorphism a of £ inherited from C which ensures Ef = Ei and Ff = Fi. 
The orbit space Ch = C/a is then, see again a Riemann surface, namely a hyperelliptic 
surface given by 



/4 



(K - iz)(K + iz) l[(K - Ei)(K - Fi) 



(2.21; 



Thus it is possible to construct components of the matrix $ on Ch which makes it possible 
to use the powerful calculus of hyperelliptic Riemann surfaces. These functions may be 
lifted to C. As we will show in the following, it is possible to construct a matrix $ on £ in 
accordance with the conditions of Theorem 2.1 by projecting onto this surface. 

3 Hyperelliptic solutions of the Ernst equation 



3.1 Theta functions asscociated with a Riemann surface and the 
Riemann— Hilbert problem 

In this section, we want to give an explicit construction of the matrix <3> in accordance with 
Theorem 2.1. Condition II can be used to construct solutions by prescribing the poles, 
essential singularities and cuts of $ which is equivalent to the solution of a generalized 
Riemann-Hilbert problem for the matrix $. The investigation of such matrix Riemann- 
Hilbert problem turns out to be rather difficult and is not yet fully done (in general it can 
be merely reduced to the solution of a linear integral equation, see e.g. [p2]). Therefore 
we will use here a different approach. The occurrence of the monodromy matrix suggests 
that it might be possible to construct a matrix $ on the Riemann surface £ of the previous 
section. The additional freedom we thus gain is used to restrict the problem to a scalar one, 
namely to a Riemann-Hilbert problem for one component of $ on the hyperelliptic surface 
Ch obtained from L by factorizing with respect to the involution a. We impose the reality 
condition Ei, Fi G K or E^ = Fi on the branch points in order to satisfy the reality condition 
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of Theorem 2.1. Then we construct the whole matrix $ in accordance with this theorem. In 
fact it was shown in [|17]] that all matrix Riemann-Hilbert problems with rational jump data 
are gauge equivalent to scalar problems on a suitably chosen hyperelliptic surface. Thus the 
limitation to the scalar case is only a comparatively weak restriction which allows, as we will 
show below, for an explicit solution of the problem in terms of theta functions. 

For the moment, we fix the physical coordinates z and z in a way that p / and that 
— iz and iz do not coincide with the singular points of $ in order to ensure that the first 
condition of Theorem 2.1 is valid. In the next section we study the dependence of the found 
solution on z and z. In order to give the solution to this special case of the generalized 
Riemann-Hilbert problem, we use the theory of theta functions associated to a Riemann 
surface (see |[23| ) and the solution of the Riemann-Hilbert problem on a Riemann surface, 
as given in p5|. As we will need only hyperelliptic Riemann surfaces of the form (12.21 ), we 
restrict ourselves to this case. 

Let us denote by (ai, . . . , a g , b\, . . . , b g ) a basis of the first (integral) homology group 
Hx(Ch) of Ch (see the picture below) where the cuts are either between real branch points 
(which are ordered E^+i < F^ + \ < . . .) or between Ei and Ei (for the moment we ignore the 
case that more than two branch points may have the same real part). 




Figure 2: The homology basis for Ch- 



dujj 



2vri5. 



Let {dcoi} denote a basis of H 1 (Ch) such that 
differentials, we define the Abel-Jacobi map of Ph € Ch by oj{Ph 
with Pq G Ch fixed. In the following, we will always choose Pq = —iz. 
We define a g x g matrix 77 - the Riemann matrix - with elements 7r 



Using these normalized 
- (fg d^, . . . du g ) 



i.i 



§ dujj . This 



matrix is symmetric and has a negative definite real part. These properties ensure that the 
theta function with integer characteristic „ defined by 



e 



(x,n) 



N 



(3.1) 
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with x G C 9 and a, (3 G Z 9 , where (•, •) denotes the Euclidean scalar product (N,x) = 
J2i=i NiXi, is an analytic function on C 9 . A characteristic is called odd if (ol,(3) 7^ mod 2. 
The Riemann vector is denoted by Kp. 

The reality condition on the branch points implies for the theta function G with charac- 
teristic q , the Riemann theta function, 

Q( x ) = 6(x + ivrA) , (3.2) 

where Aj = 1 if and Fi are real and Aj = otherwise. 

We recall that a divisor 21 on a general Riemann surface E 9 is a formal symbol 21 = 
ni Pi + • • -ntPk with Pi G S g and rij G Z. The set of divisors admits a partial ordering 
and we may associate to a meromorphic function / a divisor (/), a principal divisor. Using 
this partial ordering we define for a divisor 21 a vector space £(21) consisting of all principal 
divisors not less than 21. 

A Riemann-Hilbert problem can be stated as follows: let T be a piecewise smooth contour 
on Lh. Let A = t\ + • • • + t r be a divisor on T consisting of a finite number of pairwise 
different points subject to the following condition: T \ A decomposes into a finite set of 
connected components {Tj} (j = 1, . . . , N), each of which is homeomorphic to the interval 
(0, 1). We call Tj a curve of the contour T. Each Tj has a starting and an end point, given 
by two points of A, where the starting respectively end points may also coincide. We define 
the function a(t, Tj) on T by 

a(t '^ ) = {jo f the G rw4e ' (3 " 3) 

(j = 1, . . . , N). On each curve Tj let there be defined a Holder-continuous function Gj(t), 
which is finite and nonzero. We denote 

N 

G(t) = <*(t, T j) G i( t ) > t G T \ A. (3.4) 
3=1 

Let there be given a divisor 21 of degree m, consisting of points of the divisor A, taken 
at arbitrary degree. Let on Ch \ T be given another divisor *B of degree n. Now we can 
formulate the 

Homogeneous scalar Riemann-Hilbert problem: 

Give a function tj) with the properties 

f + (i) = G(i)f(t) , (3.5) 

with (VO g Livr 1 ^- 1 ). 

The solution of this problem was given by Zverovich [ 25) . A key point in the construction 



is the introduction of an analogue to the usual Cauchy kernel. This Cauchy analogue is given 
by a normalized (i.e. all a-periods are zero) differential of the third kind with poles at P# 
and Pq and is denoted by dup H p (r). With the Cauchy analogue at hand the solution to the 
above Riemann-Hilbert problem is given by 

4(Ph) = e* (P * } , (3.6) 
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with 



*{Ph) = ^~J InG(r) d^ Po (r) (3.7) 



where Pq ^ T, and the integration goes over all curves Tj of the contour T, and where we 

N 

have put lnG(t) = J] a(t, IT,-) lnGj(t). The 6-periods Ui of ^ are given by 

Ui = — I \nGduji . (3.8) 
27ri Jr 

Applying the Plemelj formulae to the function $?(Ph), we find 

</> + (t) = G(t)i>-(t) , t G T \ A \ |J ai . (3.9) 

i=l 

For more details on the solution to the scalar Riemann-Hilbert problem, the reader is referred 



to [||], H or |27j. 



3.2 Solutions to the Ernst equation via the scalar Riemann- 
Hilbert problem 

With the help of the results of the previous section, we are now able to state the follow- 
ing theorem, which gives the solution to the generalized Riemann-Hilbert problem on the 
hyperelliptic Riemann surface Ch- 

Theorem 3.1 Let Pq be a fixed complex constant (p ^ 0) not coinciding with the singular- 
ities of ip or the branch points Ei or F^. Let Q(Ph) be a linear combination of normalized 
Abelian integrals of the second kind (with singularities p ^ Ei and p ^ Fi, independent of z 
and z) and third kind (with in addition singularities at all real branch points with residues 
satisfying Q(Ph) = Q(Ph)- Then the solution to the generalized Riemann-Hilbert 
problem on the real hyperelliptic Riemann surface Ch is given by 



r 



, . Q(lu(Ph) -w(D) + u + b-K R ) I . 1 /• ^. NJ 
'• (/ " } = * e(J(P„)-.(D)-K R ) 6XP 1 °™ + 27i J ln G ^^(t) 

(3-10) 

where D = Pi + • • • + P g is a fixed non-special divisor on Ch which is subject to the reality 
condition: either Pi 6l or with Pi G D we have Pi G D or Pi is a branch point Ei or Fi. b 
is the vector of b-periods of Q with components 

tt , (3.11) 



i = l,...,g, the Ui are given by h3.8J , and ipo is a normalization constant. The paths of 
integration have to be the same for all integrals. 
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Proof : We want to prove that i/j(Ph) is a single valued function on Ch- If we choose a 
different path of integration for the integrals in the exponent and the map u(Ph) and denote 
the corresponding integrals by a prime, the primed and unprimed integrals are connected 
via 

Q'(P H ) = n{P H ) + j> dtt , (3.12) 
e 

(similarly for the other integrals) where S is a closed contour on Ch which may be decom- 
posed in the homology basis as follows 

9 9 

£ = ^2 mi ai + ^2riibi , (3.13) 

i=l i=l 

with 777, j, Tij G Z. Then we have, e.g. for Q and u 

n(p H ) - n(p H ) + £ ™A = n(^) + 6) , 

cj(P h ) -> w(Pir) + 27riM + 77 JV , (3.14) 

where M = (mi, . . . ,m g ), N = (ni, . . . ,n g ) G Z 9 . Under this transformation, the original 
quotient of theta functions in ( p.!0|) will be multiplied by 

exp (-<#,&)) , (3.15) 

but this term is just compensated by the contour integral over S in the exponent. The same 
argument holds for the line integral over the contour T since the Ui are its 6-periods. This 
shows that ip(Pn) is a single valued function on C H . 

^From the properties of the theta function, we also find that ip(P) has g simple poles at the 
points Pi, . . . ,P g and g simple zeros. Additional poles, zeros and essential singularities can be 
obtained by a suitable choice of Abelian integrals of the second kind (essential singularities) 
and third kind (zeros and poles). We remark that the assumption Q(P) = Q(P) had to be 
introduced in order to satisfy the reality condition of Theorem 2.1. □ 



Remark : Without loss of generality we can choose D to consist only of branch points since 
D gives the poles of ^ due to the zeros of the theta function in the denominator. This can 
always be compensated by a suitable choice of the zeros and poles of $ which arise from 
the integrals of the third kind in Q. All Pj G D shall have multiplicity 1 and be chosen in 



a way that 



(x) with 



u)(D) + lo(Pq) + K R has the same reality properties as the 



Riemann theta function Q(x) 



Our next aim is to define a matrix valued function $(P) on C, satisfying the conditions 



of theorem |2J], with the help of the above solution to the scalar Riemann-Hilbert problem 
on the hyperelliptic surface Ch- To this end we define a further function on Ch by 



x(Ph) = Xo 



G(u(P H ) +u- cu(Pq) - u(D) - K R ) 
e(u(P H )-u;(D)-K R ) 



exp 



— f \nGdujp H p 

iTXl JT 



(3.16) 
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where Xo is again a normalization constant. It can be easily seen that the analytic behaviour 
of x{Ph) is identical to that of i[)(Ph), except that it changes the sign at every ci-cut. x is thus 
not a single valued function on Ch- However, it is single valued on C which can be viewed as 
two copies of Ch cut along P , P and glued together along this cut. We define the vector X 

on C by fixing the sign in front of x m the vicinity of the points P$ = (K , 0, ±p,(K )) e C, 

X &=(±tfP))> P ~ P f- (3 - 17) 

With the help of this vector, we can construct the matrix $ on £ via 

$(P) = (X(K,fM>(K),+ii(K)),X{K,iM>{K),-ii(K))) (3.18) 

where the signs are again fixed in the vicinity of Pq\ Notice that this matrix consists of 
eigenvectors of the monodromy matrix, LX(K, fi (K),±ft(K)) = jiX(K, ^(K), ±jj,(K)) if 
L is written as L = fiQ / yQ~ l where 7 is the matrix from [2.1| . 

It may be readily checked that this ansatz is in accordance with the reduction condition 
( |2.1(J| ) (this is in fact the reason why one has to define the function x i n the way ( |3.16| )). 



The behaviour at the singularities is as required in condition II: For the contour T and the 
singularities of the Abelian integrals Q, this is obvious. At the branch points P, and Fi, one 
gets the following behaviour: at points Pi of the divisor D, the components of $ have a simple 
pole, and the determinant diverges as (K — Pi)~z, if this branch point is not a singularity 
of an integral of the third kind in Q or lies on the contour T. If the same condition holds at 
the remaining branch points, the components are regular there but the determinant vanishes 
as (K — Pi) 2. If the branch points coincide with one of the singularities of the integrals in 
the exponent in (|3.10|) , this merely changes the singular behaviour of $ and its determinant 
there. Condition II of theorem [O] is however obviously satisfied. 

Since $ in (|3.18| ) is only a function of P, it will not be regular at the cuts [Pj, Pj]. At the cl- 



outs around non-real branch points, we get $~ = $ + a 1 | a -, whereas we have &~\ ai = — & + \a i o'2 
at the a-cuts around real branch points. The logarithmic derivatives of $ with respect to z 
and z are however holomorphic at all these points. One can recognize that the behaviour 
at the non-real branch points is related to a gauge transformation of the form ( 2.12|) . This 



means that one can find a gauge transformed matrix $' that is completely regular at these 
points if the integrals in the exponent are regular there. With 

ai = i(l + A), a 2 = l -{l-\) (3.19) 



and A = U!=i \j fzf where D = £f=i P u this may be checked by direct calculation. The 
real branch points, however, cannot be related to gauge transformations. 

Normalizing ip and x (if possible) in a way that ipi 00 ^) = 1 an d x{°°h) = ~ 1) one 
can see that $ is then in accordance with all conditions of Theorem 1 since the reality 
condition follows from the reality properties of the theta functions and the Riemann-Hilbert 
problem. The fact that $ is at least differentiable with respect to z and z at points where 
Pq does not coincide with the singularities of the integrals in the exponent or the remaining 
branch points of Ch follows from the modular properties of the theta function. Let the 
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paths between [P , oo~] and [Pq, oo + ] be the same in all integrals and let them have the 
same projection into the complex plane (i.e. one is the involuted of the other). Then the 
results may be summarized in 



Theorem 3.2 Let 6 



(a>(oo ) + u) 7^ 0. Then the function 



e 


a 

A 


(co>(0O + 


) + u + b) 


e 


a 

A 


(cg>(00~ 


) + u + b) 



f(z,z) 



is a solution to the Ernst equation. 



exp \ fi(oo + ) - fi(oo ) + -^7 [ \nG(r)duj 

Z7T1 J 



T 



(3.20) 



Remarks : 



1. In the case g = the Ernst potential ( |3.20|) is real, / = e 2U . This means that U is a 
solution to the axisymmetric Laplace equation and belongs therefore to the Weyl-class. 
For g > 0, there are no real solutions other than / = 1 which describes Minkowski 
space. 

2. The multi-black- hole solutions which can be obtained via Backlund transformations 
(see e.g. H28fl ) are contained in the class ( |3.20|) as the limiting case that the branch 
points Ei and Fi coincide pairwise. In this limit, all branch points become double 
points and the theta functions break down to purely algebraic functions. Notice that 
the analysis of / at the branch points in the following section always assumes a regular 
surface. The obtained results for the regularity of / do not hold in this limit. 

The above explicit construction of the solutions makes it possible to derive useful formulae for 
the metric function a and the derivatives of the Ernst potential. Let Jp H+S dui = g^d + o(5) 
where 5 is the local parameter in the vicinity of Ph 6 C-h- We define the derivative 



z 

i=l 



D Ph Q(x) = ^2gid Xi e(x). 
Using (EuaD and (|3~T0|) , (|3~T6|) , we get 



(3.2L 



(a — ao )e 



2U 



iDoo- In 



e 


a 

A 




duo + u 


+ 6) 


e 


a 

A 




du + u 


+ 6) 



(3.22) 



^From the linear system ( |2.5 a|) and ( [2.5 bp , we obtain with ( |3.10|) and ( |3.16| ) 

fz 

f + f 



i e 


a 

A 


(u + b-uo(p ))e 


a 

A 


(m + 6 + u;(oo )) 


2^P -P Q 


a 

A 


(u + b)Q 


a 

A 


(u + b + Lu(oo~) -w(P )) 



X 



D Po In 9 



(u - u(P )) + I Po 



(3.23) 
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and 



f* 



f + f 



i e 


a 

A 


(u + b)& 


a 

.P. 


(u + 


M 


- u(0O ) - Uj(P )) 




a 

A 


(u + b-u(p ))e 


a 

A 


(u + b + u;(oo - )) 



D Po In 6 



x 



(u-u{P )) 



(3.24) 



where Ip is the linear term of the expansion of the integrals in the exponent of ( 3.10|) in the 
local parameter around Pq. 



3.3 Finite gap solutions and Picard— Fuchs equations 



The original finite gap solutions of ]14| are those among ( |3.20| ) without the contour integral 
(in our notation only an arbitrary linear combination of Abelian integrals of the second and 
third kind fl). They just correspond to the so called Baker-Akhiezer function (see ||13|| ) for 
the Ernst system. This function that has essential singularities and poles gives the periodic 
or quasiperiodic solutions to the integrable nonlinear evolution equations. There the essential 
singularity is uniquely determined by the structure of the differential equation. In contrast 
to these equations, the solutions ( |3.2(J| ) are in general neither periodic nor quasiperiodic, 
and the essential singularity can be nearly arbitrarily chosen. The form of the solution to 
the Riemann-Hilbert problem shows that one might even think of "putting the singularities 
densely on a line and integrate over the integrals with some measure": an Abelian integral 
Q p of the second kind with a pole of first order at p can be used as an analogue to the 
Cauchy kernel. A contour integral over this kernel with some measure, J r lnGQ p dp, is thus 
just another way to write down the solution to a Riemann-Hilbert problem on a Riemann 
surface. 

In studying the boundary value problem for the rigidly rotating disk of dust, Meinel and 
Neugebauer [HJ observed that it is possible to obtain solutions to the Ernst equations via 



/ = exp(t / 

\m=l JE 



c m K9dK 



(3.25) 



where the divisor C = J2m=i C m is determined by 

r c m K { dK 



E 

m=l 



(3.26) 



0,1, 



,g— 1), i.e. as the solution of a Jacobi inversion problem. The Jj are (in the 
absence of real branch points) real solutions to the axisymmetric Laplace equation which 
satisfy the recursive condition, 



il 



n+l,z 



zl, 



1 



(3.27) 



The relation to the class obtained in theorem |3^ is the following: The integral of the 
third kind in ( |3.25 ) can be expressed by the help of a formula in |33| via theta functions. 
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Equation (|3.26|) ensures that the resulting expression is independent of the chosen integration 
path which is shown in the proof of the theorem. Thus the (obtained from the Jj by 
normalization) are as in our case the 6-periods of the integral I g in the exponent. In fact, it 
was shown in [[U| that one of these periods, say Ii, can be chosen as an arbitrary solution to 
the axisymmetric Laplace equation. The other periods as well as the integral in the exponent 
then follow from differential identities plus boundary conditions. 

The underlying reason for this fact is that the Ernst potential / is studied on a family of 
Riemann surfaces parametrized by the moving branch points — iz and iz. The periods on this 
surface (i.e. integrals along closed curves) are subject to differential identities, the so called 
Picard-Fuchs equations. It is a general feature of the periods of rational functions p9|, [JT 



that they satisfy a differential system of finite order with Fuchsian singularities. An elegant 
way to find the Picard-Fuchs system explicitly is via the notion of the Manin connection 
in the bundle H^ R (T> g ) — > E fl , see f32|. The investigation turns out to be particularly 
simple if one uses the following standard form of the (hyperelliptic) Riemann surface E s (all 
hyperelliptic surfaces of genus g are conformally equivalent to this standard form) 



23 



[x 



-z)H(x-E i ) = (x-z)P(x) 



(X 



29 

j=0 



(3.28) 



where the E { do not depend on z. Using j = dx/y, j\ = xj , . . . , j2 9 -i = % 29 1 jo as the basis 
for the de Rham cohomology H^(E g ) we obtain for the matrix M™ (m, n — 0, . . . , 2g — 1) 



of the Manin connection (defined by 



din 



lvl n Jmj 



Ml 




m + l)a m+ i + z 



-m—l 



a i z * 

j=0 



for < m < n, 



2g— 1— m 



(3.29) 



m + l)a m+ i - a-m+i+jZ 3 J for n < m < 2g - 1 

j=0 



One finds that the periods satisfy a similar recursive condition as ( |3.27 ). An analogous 
consideration can be performed for the ^-dependence of the periods. One finds that the inte- 
grability condition of the Picard-Fuchs systems is just the axisymmetric Laplace equation. 

On the other hand, with the help of some boundary conditions (for instance at \z\ — > 00), 
the I n can be uniquely determined from the above system (|3.27p . Thus the class of solutions 
discussed by Meinel and Neugebauer may be phrased in the following form: if an arbitrary 
solution of the Laplace equation is given, one can calculate the functions I n with ( |3.27| ) 
and the boundary condition, and ends up with a solution to the Ernst equation of the form 



4 The singular structure of the Ernst potential 



The construction of the solutions in the previous sections with the help of Theorem |2.1| also 
indicates where the resulting Ernst potential ( |3.20| ) may be singular: only at points P = ~ lz 
where the conditions of Theorem 12.11 do not hold. Notice that these conditions are sufficient 



for the regularity of / at all other points. It may turn out though that the Ernst potential is 



15 



perfectly regular at points where Theorem 2.1 is not fulfilled, e.g. in the case of singularities 



that are pure gauge. We will therefore discuss all possible singular points of the solutions 



It is very helpful that this whole discussion can be performed on the Riemann surface 
Ch where one can use the powerful calculus on hyperelliptic surfaces. One does not have to 
work on the four sheeted surface C whose introduction was necessary for the construction 
of the solutions, and which provides an understanding of the mathematical properties of the 
Ernst equation. Since we will work from now on on Ch only, unless otherwise noted, we will 
drop the index H at points P G Ch- 

The possible singularities of / can be directly inferred from the potential in the form 



( 3.2C ). The Ernst potential will be singular at the zeros of the denominator. It is possibly 
not regular at the points where Pq is identical to the singularities of Q or is on T. Critical 
points of a different kind are the branch points Ei and Fj. If P coincides with these points, 
the Riemann surface Ch becomes singular. Something similar happens at the axis where 
the branch points Pq and Pq coincide. This is a reminiscent of the singular behaviour of 
the three-dimensional Laplace operator on the axis in the axisymmetric case. The main 
aim of the following analysis is to single out a class of solutions that may be interesting 
in the context of boundary value problems for the Ernst equation that describe e.g. the 
exterior of a body of revolution. Thus we will not study the nature of the singularities (e.g. 
curvature singularities) but single out a large class of solutions where the Ernst potential is 
only discontinuous at a (closed) contour that could be identified with the surface of a body. 



Zeros of the denominator Zeros of the denominator of ( |3.20|) will lead to singularities 



in the spacetime. From condition IV of Theorem 2.1 it follows that these are just the 
points at which the matrix $ cannot be normalized in the required way. This leads to the 
transcendental condition 

9 a (u(oo~) +u + b) ^0, (4.1) 
_P _ 

if one wants to exclude these zeros of the denominator. We will show in the next section 
how the points at which 6 ^ (co>(oo~) + u + b) = can be found as the solution of a set of 
algebraic equations. 

Essential singularities The integrals of the third kind occuring in Q are nothing but a 
particular case of line integrals over contours with constant jump function G(t). Therefore, 
we are left with the investigation of the integrals of the second kind at this point. Since the 
theta functions in ( 3.20|) are regular as long as the Riemann surface Ch is, we are left with 



the exponent if (f4.1|) holds. For the behaviour of the exponent, we get the following 



Proposition 4.1 The Ernst potential ( \S.2Q ) has an essential singularity at the points where 



Pq coincides with the singularities of the integrals of the second kind on Ch- 

Proof : The exponent has by construction an algebraic pole there and, therefore, the Ernst 
potential has an essential singularity. □ 
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An essential singularity of the real part of the Ernst potential corresponds to a line singularity 
of the metric function U. In the context of exterior solutions for bodies of revolution we 
are interested in, there seems to be no situation where such a line singularity in a spacetime 
might be interesting. 



Contours In the case that Pq lies on a contour Tj but not on an endpoint of T,, on the 
axis, or on one of the branch points Ei or F{, it can be easily seen that the integral in the 
exponent as well as the 6-periods U{ are bounded since G is Holder-continuous, finite and 
non-zero on T. At the endpoints of the contours singularities may occur. The value of 
these integrals at the remaining points will however not be the same in general if the contour 
is approached from one or the other side. This can be seen from the following fact: The 
point Pq is a branch point of Ch- If it lies on the contour T, care has to be taken of the 
sign of the root whilst evaluating the integrals of the form J n = J r lnG(r) 1 -^ 2: . The decisive 
factor is $ = (K - P )(K - P ). We get for K e V with K = K x + iK 2 , and K U K 2 eR 
for the imaginary part of Hq, 

c^ = ±-^sgn ((Jfc - ()K 2 ) >/|a*8|-R(m8), (4.2) 

i.e. the sign of the imaginary part of p depends on the sign of K\ — C. Thus the value of 
the integrals will in general not be the same whether the contour is approached from the 
interior or the exterior region. This reasoning does not work for points K = K not on the 
axis (K 2 7^ 0) with K\ = 0. There the imaginary part of pfj is zero which means that po is 
either purely imaginary or real in the vicinity of Pq = Kq depending on the sign of K 2 — (■ 
We conclude that the integrals over T with Pq e V have the form J = J 1 + sgn(e) J 2 (where 
the J 1 are independent of e which indicates if the contour is approached from the interior 
or the exterior) which implies that the limiting value of the Ernst potential calculated via 
( p.20| ) exists but depends on e. Therefore, we have proven the following 

Proposition 4.2 Let Pq lie on the contour V but not on the axis, on one of the branch points 
Ei or Fi, or at an endpoint ofTi. Then f will, in general, have a jump at T. The limiting 
value of f will exist and be Holder-continuous there, f may be singular at the endpoints of 
the iy 

Thus the Ernst potential will be finite but discontinuous at a contour T z in the (p, £)-plane 
given by Pq G T which means that the solution to the vacuum equations will not be regular 
at a surface in the (p, (, 0)-space. If this surface is closed, it can possibly be identified 
with the surface of a body of revolution. The interior of the body is supposed to be filled 
with matter. Therefore the vacuum solution is only considered in the exterior (that contains 
z = 00); it is not regular at the boundary to the matter region. 



The axis The axis is a double point on the Riemann surface Ch since two branch points 
coincide. In this case, all quantities may be considered on the Riemann surface E' given 
by ji' 2 = nf=i( r — Ei)(r — F^, as was shown by Fay |34 



Let a prime denote here and in 
the following that the primed quantity is taken on This surface is obtained from Ch 
by removing the cut Pq, Pq . For the analysis of the axis, we will use a slightly different 
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cut-system than the one introduced in the previous section: we take a closed curve around 



unchan 



in the H — sheet as the cut a g . All 6-cuts shall begin at the cut [Ex,Fi\. The rest is 
ged. This implies for the characteristic of the theta function that it has the form 



a' 1 

P' e 



(4.3) 



where e — 0, 1 and a' { = 0. 

Since the expansions of all characteristic quantities of the Riemann surface are smooth 
in p except ir gg which is divergent as In p for p — > 0, it follows that the Ernst potential has a 
regular expansion in p. For points Pq not coinciding with real branch points or singularities 
of the exponent in (|3.20|) , the Ernst potential is thus at least C 3 . It follows from a theorem of 
Miiller zum Hagen |33] that it is therefore analytic. Consequently it is sufficient to calculate 
the limiting case. If this is well defined, the Ernst potential is regular at these points of the 
axis. The differentials of the first kind for p = turn out to be 



duJi = da;', 



,g-l, du g 



-dc^ +c _ 



(4.4) 



where du'^+f- is the normalized differential of the third kind on £' with poles in ( + and ( 
This implies for the 6-periods 



71 



n, 



C+ , 
2 In p + reg. terms 



l,...,g- 1, 
i = 1, ■ ■ ■ ,9 ~ 1, 



(4.5) 

(4.6) 
(4.7) 



Since ir gg diverges, the theta function will break down to a sum of two theta series on E' (in 
the case of genus g — 1, the surface E' has genus 0; the formula below can however be used if 
one replaces the theta function G' simply by a factor 1 which means that the axis potential 
can be expressed in terms of elementary functions in this case). We introduce integrals u)'(P) 
of the first kind with the property u'(Ei) = 0. The differential dc^oo+oo- on the axis becomes 
du' + . In the case of the contour integrals one has to observe that an additional factor 



sgn(i^! — £) in the notation of (|4.2|) occurs for the same reasons as there. Since the Abelian 
integrals of the second kind can be obtained from the integrals of the third kind by a limiting 
procedure, the same holds for these integrals and their 6-periods. With the above settings, 
we obtain for ( |3.20| ) 



./ 



0' 


V" 

.P. 


(</ 


~ + + + u' + V) 


+ (_l) £e xp(-K(oo+)+ % + 6 fl ))e' 


a'~ 

.p. 






oo + + u > + b>) 


6' 


'a 1 ' 

3'. 


(" 


™ + + - u< - b>) 


+ (-1) £ exp(-(^(oo+) u g b g ))& 


'a 1 ' 

.P. 




p - u' - b') 



exp + 



J In G(<r)da4 +00 - (r) + b g + u g J 



(4.8) 



It can be seen from the above formula that the limiting value of / exists even if u g diverges, 
provided (|4.1| ) holds (/ will be Holder-continuous if u g diverges). The Ernst potential will 
however have an essential singularity at the real singularities of Q. We can summarize the 
above results. 
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Proposition 4.3 Let condition ( \4-A ) hold. Then the Ernst potential is regular on the axis 
except at the points where P coincides with singularities ofQ, points ofT , and branch points 
E F 



Remark : Though / is Holder-continuous even if u g diverges, it is interesting to note for 
the following when this will be the case. Obviously this can only happen at the real points 
of r. It can be seen however that u g is always bounded at these points due to the reality 
condition, unless they are endpoints of Tj (this would lead to a conic singularity on the axis). 



Real branch points If P coincides with a real branch point E^ or F i} this will be a triple 
point on C H . We get the following. 

Proposition 4.4 At points where Pq coincides with the real branch point E g , the limiting 
value of f exists. The Ernst potential is in general not differentiable there. 

Proof : We use the same cut system and the same notation as on the axis. Put Pq = E g + x 
with x = 5e 1( ^ and (ft G R, 5 G M + . In order to expand / in powers of x and x, one has to 
consider the a-periods, in particular 

dr 4 d dt 



* 3 V \jx(F g - E g ) •>! yj(l - t 2 )(l - kH 2 )(F g -E g + xt 2 )^"(F g + xt 2 ) 

■ 4 = : M (k) + 0(5)) (4.9) 
x(F g - E g )^»(F g ) 



where k = e lfli , where K(k) = K(yl — k 2 ) and K(k) are the complete elliptic integrals of 
the first kind, and where // /2 (t) = nf=i( r — E,j){r — Fj). It can be seen from ( |4.9|) that the 
a-period has an expansion in powers of y/5. The coefficients of the expansion in y/x, and 
y/2 are ^-dependent, since the modul of the elliptic integrals is just k = e 1 ^. This implies 
for the differentials of the first kind duji = du^ + 0(VS) for i = 1, g— 1, and du g = du; 9 _ 1 . 
Similarly duj^+oc- = duj' 00+Q0 - + O (yty ■ We get for the 6-periods, 

^ = - 2T H( 1 + (^)) • (410) 

whereas ^( g -i) g = 0(\f8) and ixij = tt'^ for i, j = 1, . . . ,g — 1 in the limit. Thus / can be 
expanded in yfx and Even in case that only integer powers in the expansion occur, the 
coefficients will be in general (^-dependent. Though the limiting value of / at Pq — E g exists, 
/ will in general not be differentiable at this point. □ 

This implies that the real branch points are singular points on the axis, possibly topological 
defects in the spacetime, see [frjj]. They should not occur in the context of exterior solutions 



for bodies of revolution we are interested in. 
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Non-real branch points If P coincides with a branch point E g = F g , the points E g and 
F g will be double points on Ch- Thus the situation is similar to the one on the axis with the 
only exception that one ends up here with two double points. As on the axis, it is convenient 
to consider all quantities on a Riemann surface S" given by /i" 2 = nf=i (K ~ Ei)(K — Fj) 
where the double points are removed. All quantities with two primes are understood to be 
taken on this surface. We use the following cut system: let a g _i be the circle around [P , E g ], 
and a g the circle around Po,F g , both in the plus sheet. The remaining cuts are as on the 
axis, i.e. all 5-cuts start at [Ei,Fi\. As on the axis, we get 



Proposition 4.5 Let hold, and let E g = F g ^ F . Then the Ernst potential is regular 
at the point P = E g . For E g 6 T , / is in general Holder-continuous at P = E g . 

The proof is similar to the one on the axis and basically uses again results of Fay MM . 

Proof : The case g = 1 may be checked directly with the help of the standard theory of 
elliptic theta functions (see e.g. |3(|). For g > 1 with the cut system in use and Pq = E g + x, 
where x is chosen as in the case of the real branch points, the differentials of the first kind 
have a smooth expansion in x and x. In contrast to the case of real branch points, the 
coefficients in the expansion are (^-independent. The differentials dcUj become in leading 
order the differentials of the first kind du" on £". The differential du; g _i becomes in the 
limit the differential — du/' + „_, and similar for du) at F a . The differential of the third 

hj g h g y y 

kind becomes d^oo+oo- = da;^ +no _. All these differentials have coefficients in the x and 
x expansion that cor 
as may be checked b 
i,3 = !, --,9- 2 and 



x expansion that contain Abelian integrals of the second kind with poles in E g and F^ 
as may be checked by direct calculation. This implies for the 6-periods that 7Tjj = tc'-j for 



TRs-iKs-i) = ir gg = 2\n5 + 
TTi^x) = -2u"{E+), 

7t ig = -2^"(F 9 + ), (4.11) 

whereas 7c^ g _i^ g is finite in the limit 5^0. If E g ^ T, the Ui as well as the Cauchy integral 
in the exponent have a smooth expansion in x and x with finite coefficients. The theorem 



of |35j then guarantees regularity if the limiting value that may be calculated as on the axis 
exists. The theta function on Ch breaks down to a sum of four theta functions on S" times a 
multiplicative factor. If E g G T, the coefficients in the expansion of / in x and x will diverge 
which implies that / is possibly not different iable there though the limiting value exists if 
m) holds. □ 



5 Metric functions and ergospheres 

In the previous sections we have made extensive use of the complete integrability of the 
Ernst equation to construct a large class of solutions. To discuss physical features of the 
resulting spacetimes however, it would be helpful to have expressions in closed form not only 
for the Ernst potential but for the metric functions, at least for the functions e 2U and a that 
can be expressed invariantly via the Killing vectors. It is a remarkable fact already noticed 
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by Korotkin [13 that the metric function a can be related to derivatives of the matrix $ 
without solving the differential equation ( [2.2|) . In the following we will show that a theta 
identity of Fay []33| can be used to go one step further to obtain a formula for a that is free 
of derivatives. The same identity leads to a simplified expression for the metric function e 2U 
that can be directly used to identify ergospheres in the spacetime. 

Fay's trisecant identity establishes a relation between four points Ai, . . . , A± on a Riemann 
surface, in our case in arbitrary position (see e.g. | ]37[| , |38|| ). Let x be an arbitrary g- 
dimensional vector. Then the following identity holds, 



0(x)0 ( x + 1^ 3 du + du \ - exp (^AiaJaI) 6 f £ + " did J 6 f x + < L: 



A-2 



A-2. 



- exp {SIa 2 aAm) 9 + f A 3 duJ ^j + J A 4 dcJ ) = 



(5.12) 



where e.g. ^AjaJ^o denotes the integral of a normalized differential of the third kind with 
simple poles at Ai respectively A 4 with residues +1 respectively —1 along a path from A 2 
to A 3 . For a geometric interpretation of this identity in terms of the Kummer variety see 
for an interpretation via generalized cross ratio functions see 24]. The strength 
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of the above identity arises from the fact that it holds for points A^ in general position. By 
a suitable choice of these points, we obtain for the metric function e 2U , the real part of the 
Ernst potential, 
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(u + b + lo(oo-) + lo(P ))Q 
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(u + b + u(oo )) 



(5.13) 



where / denotes the integral in the exponent of ([3.20 ). This formula makes it possible 
to identify directly the zeros of e 2U which give the ergospheres, the limiting surfaces of 
stationarity (inside these surfaces there can be no observer at rest with respect to spatial 
infinity). Since the exponent of the integral of the third kind in ( p.!3[ ) in front of the fraction 
cannot vanish, the necessary condition for ergospheres is 







{u + 6)0 



{u + b + uo{P Q )) =0. 



Defining the divisor A as the solution of the Jacobi inversion problem 

uj{A)-uj{D) = u + b, 



(5.14) 



(5.15) 



we find that an ergosphere can occur if Pq or Pq are in A. It is however possible that the 
denominator of ( |5.13| ) vanishes at the same time which would imply a violation of ( |4.1| ) (and 
thus a singularity of the spacetime). Summing up we get 



Proposition 5.1 /. Let Pq or Pq and oo~ be in A for some Pq, then condition ( \j.l\ ) is 
violated and the Ernst potential is singular at these points. 

II. Let P or P but not oo~ be in A for some Pq, then the real part of the Ernst potential 
vanishes at these points which describe an ergosphere. 
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The same formula can be used on the axis where we obtain for the metric function e 2U 
in the notation of 



2U 



exp [Q 



0' 



K + V) 



0/ 



l2 



P' 



[00' 



+ u' + V) - exp(2(c; g (oo-) + u g + o 9 ))0' 



0' 



[00' 



The condition for an ergosphere to hit the axis is then 



0' 



a' 



(vf + 6') = 



+ v! + V) 

(5.16) 



(5.17) 



since the integral of the third kind in the exponent in front of the fraction cannot diverge for 
finite values of (. The interesting feature of this relation is that it is completely independent 
of the physical coordinates. This implies that if an ergosphere extends to the axis, this will 
be only possible if the metric function e 2U vanishes on the whole axis. In the case of the 
Kerr solution, the ergosphere touches the axis at the horizon. An interpretation of the fact 
that the whole axis would be singular in the present case is given in the next section where 
the above case is related to the ultrarelativistic limit in which the source of the gravitational 
field becomes so strong that it vanishes behind the horizon of the extreme Kerr metric. 

The metric function a can be calculated from ( |2.13| ) if one uses the trisecant identity 
( |5.12| ) in the limit that two points coincide. Using the trisecant identity several times, we 
get 
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(5.18) 



The constant ao can be obtained in a similar way from the condition that a = on the 
regular part of the axis (we assume here that the singularities in the exponent of (|3.20|) are 
situated in a compact region of the (p, £)-plane). Care has to be taken in the above formula 
that some of the terms in brackets explode as 1/p in the limit p — > 0. We get 
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(5.19) 



It can be seen from this formula that ao does not vanish if there are no singularities in 
the exponent (/ = u — b = 0) in which case / = 1 which describes Minkowski spacetime. 
This reflects, as already noted, the fact that a is a gauge dependent quantity. The metric 
function a however is gauge independent. In the above example of Minkowski spacetime, it 
will of course vanish in the used asymptotically non-rotating coordinates. 
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6 Asymptotic behaviour and equatorial symmetry 



Since we are mainly interested in solutions to the Ernst equation that could describe the 
gravitational field outside a compact matter source, we will study the asymptotic behaviour 
(near spatial infinity) of the solutions ( |3.20|) . It is generally believed that the Ernst potentials 
of the corresponding spacetimes are regular except at the contour T z in the (p, C)-plane which 
corresponds to the surface of the body, asymptotically flat and equatorially symmetric. 
We will investigate in the following whether it is possible to identify solutions with these 
properties in the class ( |3.20| ). 



Asymptotic behaviour Asymptotic flatness implies that the Ernst potential is of the 
form f — 1 — 2m/\z\ + o(l/|z|) for \z\ — > oo where m is a positive real constant. A complex 
m is related to a so called NUT-parameter that is comparable to a magnetic monopole. 

The asymptotic properties of the solutions (|3.20|) can be read off at the axis. Notice that 
the d.ui[ are independent of (. For du g , we get 




du g = du4 +00 _ 1 - - Y,(Ei 7 dw Ui + o (V0 (6-20) 



where du;^ 0+ 1 is the differential of the second kind with a pole of second order at oo + . 
Furthermore it can be seen that exp(— u g (oo + )) is proportional to l/( for £ — > oo. Thus we 
get 

Proposition 6.1 Let lim T _ >00 rlnG ! (r) = on all contours that go through oo + or oo~ and 

let Q' a, (vf + b') ^ 0. Then f has the form f = 1 — 2m/ X for £ — ► oo where m is a complex 
constant. 



The proof of this proposition follows from (6.20) and (|4.8|) 



Equatorial symmetry The fact that the mass is in general complex implies that the class 
we are considering here is too large if one wants to study only solutions that are asymp- 
totically flat in the strong sense (m real). There is the belief that stationary axisymmetric 
spactimes describing isolated bodies in thermodynamical equilibrium are equatorially sym- 
metric. This implies for the Ernst potential /(— C) — /(C)- Solutions with this property 
always have a real mass due to the symmetry. It is therefore of special interest to single out 
equatorially symmetric solutions among those in ( |3.20|) . We get 



Theorem 6.2 Let Ch be a hyperelliptic surface of the form (\2.21\) with even genus g = 2s 



and the property p(—K, — () = p(K, (). Let T be a piecewise smooth contour on Ch such that 
with P = (K, p(K)) G T also P G T and (-K, p(K)) G T. Let there be given a finite nonzero 
function G on Y subject to G(P) = G(P) = G((—K,p(K))). If (p, p(p)) is a singularity 
of Q, the same should hold for (—p,fi(—p)). Choose a cut sytem in a way that the cuts a} 
(i — 1, . . . , s) encircle [—Fi, —Ei] and af encircle [Ei, Fj\ in the +-sheet (in the case of real 
branch points, the points are ordered in the way E^ < F^ < E i+ i < . . .; points with the same 
real part are ordered in the way 5s (Ei) < 5s(Ei) < 5s(E i+ i) < . . . which implies that Ei ^ E 
in this special case). 
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Then f is equatorially symmetric if the characteristics in the i-th position (any combination 
of the two cases is allowed) have the form 



(6.21) 



"o o" 




"o o" 


1 1 


> 






Proof : The property fj,((, K) = yu(— C, — K) on Ch makes it possible to express quantities 
on a surface with ( = —Co in terms of the corresponding quantities on the surface with 
C = Co- We have aj(—Q = ™|(C) an d b}(— £) = rfef(C) where C and — ( denote the surface 
on which the quantity is considered, and where r is the anti-holomorphic involution on Ch- 
Together with the symmetry properties of the Abelian integrals in the exponent of (|3.20|) , 
this implies that the transformation C — > — C acts as the complex conjugation together with 
a change of the upper index. Thus we have for the characteristics in ( |6.21|) /(— Q = /(C)- 

□ 



Remark : If the theta function contains only blocks of the form 





1 1 



the resulting / is just 



the complex conjugate of the Ernst potential built with the Riemann theta function. This 
means that the two cases in fl6.21|) are related through complex conjugation. It is however 
possible to combine any number of these two blocks in which case the Ernst potential cannot 
be simply reduced to the case of the Riemann theta function or its complex conjugate. In 
the case of a rotating body, complex conjugation of the Ernst potential only implies that the 
angular velocity of the body changes its sign. 

The above results suggest that it is possible to identify a whole subclass of solutions 
among ( 3.20Q that are asymptotically flat, regular except at a closed contour and equatorially 
symmetric, i.e. solutions that might describe the exterior of a rotating body and might be 
helpful in the construction of solutions to boundary value problems for the Ernst equation. 
We get 

Theorem 6.3 Let Ch be a regular hyperelliptic surface of even genus g = 2s of the form 
f$-2\) without real branch points. Let T be a closed, smooth contour on Ch such that with 
P = (K,fi(K)) E T also P = (K,n(K)) E T and {-K^n{K)) E V and E { £ V. Let there be 
given a finite nonzero function G on T subject to G(P) = G(P) = G{(—K,n{K)). Choose 



the characteristic 

theorem 6.2. 
Then 



such that it consists of blocks of the form 
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-exp < — / lnG(r)dw 

_ 7/.1 2.1X1 J 



T 



as in 



(6.22) 



■is 



1. a regular solution to the Ernst equation for P if condition (JA ) holds. 

2. in general discontinuous at T z given by Pq ET, 

3. asymptotically (\z\ — > oo) given by f = 1 — 2m/\z\ where m is a finite real constant if 



9' 



u 



4- equatorially symmetric. 
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Proof : From ( |3.20|) it can be seen that / is a solution to the Ernst equation. The regularity 
properties follow from the previous section. Asymptotic behaviour and equatorial symmetry 
follow from above. □ 



Remark : The choice of this class is mainly due to regularity requirements. If all singulari- 
ties like real branch points or the singularities of the Abelian integrals of the second kind lie 
within the contour T where the solution is not considered since the region is assumed to be 
filled with matter, they would not affect the vacuum region. However this would not enlarge 
the degrees of freedom (one real-valued function and a set of complex parameters) if one 
wants to solve boundary value problems. 

We will discuss the common properties of the solutions in this subclass in the following. 



Mass and Angular Momentum The asymptotic behaviour of the Ernst potential, / = 
1 — 2m/\z\ — 2ij/\z\ 2 + • ■ •, follows already from the axis potential ( |4.8| ). The equatorial 
symmetry implies that m and j are real constants. In the following it is convenient to 
introduce rescaled coordinates z = z/R where R is the radius of the smallest sphere that 
totally contains the contour T z , and the dimensionless quantities M = m/R, J = j/R 2 . 
This implies that the contour shrinks to a point in the limit R — > 0. If we use the differential 
operator we get for the ADM-mass with 
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Doo + 
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a 1 ' 
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2pi Jr 



]nGdu[ 



(6.23) 



A solution is of course only physically acceptable if the mass is positive. Similarly one can 
see that the angular momentum is proportional to I/O 2 



P' 



Minkowskian limit It is possible to parametrize a solution by the mass and the angular 
momentum. For M << 1 and J « 1, the solution is nearly Minkowskian. It can be 
directly seen that this limit is obtained for G — > 1. This implies that the solutions from 
above for \G\ ~ 1 are in the regime of small gravitational fields. This is also the regime of 
the Newtonian limit if the solution has one. 



Ultrarelativistic limit 

momentum diverge if O' 



P' 



It can be seen from Q6.23 ) that both the mass and the angular 
u') = but that M 2 / J remains finite in this case. This suggests 



that this divergence is best understood as the limit R — > as was already done in |39j for the 
case of the rigidly rotating dust disk: in this limit, the gravitational fields at the surface T z 
of the matter source become so strong that it is hidden behind a horizon, i.e. its coordinate 
radius tends to zero. It is suggestive to consider this limit as the ultrarelativistic limit of 
the solution. For finite R, the Ernst potential is purely imaginary on the axis which means 
that the solution is no longer asymptotically flat. If one takes the limit R — > 0, one ends 
up with the Kerr solution in this case which gives further support to the interpretation of 
this limit as the ultrarelativistic limit (this interpretation is of course only consistent if the 
limit gives the extreme Kerr solution for which the horizon is given in the used coordinates 
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by p = C = 0). If the limit R — > is taken for p = ( = with p/po and C/Po finite (this 
corresponds to an observer on the contour that vanishes behind the horizon), the resulting 
solution will not be asymptotically flat. For an observer on T z , the exterior region is in 
infinite geodesic distance, and thus completely decouples from the exterior. In the context 
of a boundary value problem this limit corresponds to a stability limit for the solution: if 
the boundary data reach a certain critical value, the solution will no longer be regular. 

Thus solutions in ( p.22| ) should be physically interesting in the parameter range < 
M/R < oo. Notice that the solutions have an analytic continuation beyond the upper limit. 



a' 

have a 



It can be seen however from ( |6.23|) that the mass changes its sign in this case since 
has zeros of first order. Consequently these 'overextreme' solutions, that do not 
Newtonian limit, are probably not physically interesting, at least they are unacceptable in 
the region with negative mass. 



7 Reduction of the Ernst potential 

The explicit form of the solutions ( |6.22| ) in terms of theta functions has a number of advan- 
tages in contrast to the linear integral equations to which the solution of boundary or initial 
value problems can be reduced in the case of integrable non-linear evolution equations^} it is 
possible to identify physically interesting features as ergospheres or the ultrarelativistic limit 
explicitly. Since the theta functions are transcendental, final results normally can only be 
obtained numerically. It is however possible to address most features like the condition for 
the occurrence of ergospheres (|5.14j) directly without having to determine the Ernst potential 



numerically in the whole spacetime. 

The numerical treatment of theta functions is comparatively simple since the exponential 
series converges rapidly due to the factor exp {^Tiij-riin^j where 3?(7Ty) < 0. It is however 
obvious that the numerics become more and more tedious the larger the genus g of the 
Riemann surface Ch is. Therefore it is an important question whether the Riemann surface 
can be reduced in physically interesting cases to surfaces of lower genus. Loosely speaking 
this is possible if there exists a special relation between the branch points (see Weierstrass' 
discussion of the case g = 2 which is referred to in |TjJ and references given therein). Since 
the branch points Pq, Pq are parametrized by the physical coordinates and can thus take 
on arbitrary complex values, such a reduction will only be possible at special points of the 
spacetime which will in general not be of special physical interest. A general reduction of 
the Riemann surface is possible if there exist non-trivial automorphisms on the surface. For 
the class of equatorially symmetric solutions discussed here, this is the case in the equatorial 
plane and on the axis. There the surfaces Lh and £' have defining equations p(K) and p'{K) 
which both depend only on K 2 . Thus on both surfaces there is the involution T defined by 
(K,p(K))^(-K,p(-K)). 

For the sake of simplicity, we will only discuss the characteristic on = f3{ = and the case 
El = —Ef (the general case can be inferred from the resulting relations without problems). 
We will concentrate on disks of radius po since they are an interesting model for galaxies. The 



*There are solutions in terms of theta functions for these equations, too. But as we have pointed out 
already, these solutions are always periodic or quasiperiodic. The solutions to the Ernst equation discussed 
here are however gauge equivalent to solutions to a linear integral equation as was shown in Jl7|]. 
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Ernst potential simplifies in the equatorially symmetric case at the disk where the boundary 
data are prescribed, what makes disks the most promising objects in the search for solutions 
to boundary value problems to the Ernst equation in closed form. We recall that the first 
solution of such a problem was found for the rigidly rotating dust disk |15| . 

In the equatorial plane (( = 0), the surface Ch is then given by p 2 (K) = (K 2 + 
P 2 ) nf=i {K 2 - E 2 ) (K 2 -El). We cut the surface as before which implies Taj = a 2 , Tbj = b 2 
and dui}{TP) = —duif(P) with P e C H . The Riemann surface Si = Ch/T of genus s is 
then given by 

rf(x) = x(x + p 2 ) f[(x - E 2 )(x - Ef) . (7.24) 

i=l 

The holomorphic differentials dvi in Ei dual to (a^, 6j) (the projection of the cuts on Ch onto 
Ex) follow from dvi = duo] — duo 2 . The so called Prym differentials dwi which change the sign 
under T are given by dw{ = duj + du 2 . They are holomorphic differentials on the Riemann 
surface E 2 of genus s with 

tily) = (y + p 2 ) fl(y - $)(v - e 2 ) , (7.25) 

which implies that the Prym variety is a Jacobi variety in this case. The Riemann matrix 
on Ch has the form 

n= 2l v n 2 -n 1 n 1 + n 2 J ' (7 - 26) 

where the II 1 are the Riemann matrices on E J respectively. The theta function on Ch thus 
factorizes into products of theta functions on the Ej, 



9(^2,11) = £e 

8 



(xi + x 2 ;2n 2 )6 



(x 1 ~x 2 ;2U 1 ) , (7.27) 



where each component of the s-dimensional vector S takes the values 0, 1. Thus the theta 
function on the surface of genus 2s can be expressed via theta functions on surfaces of genus 
s. 

In the case of the Ernst potential ( p.22| ), further simplifications follow from the fact that 



oo is a branch point of E 2 . For the contour integrals u, we obtain for disks 

u v = — / InGdv, u w = sgn( — / InGdw , (7.28) 

7T1 JF V 7T1 Jt w 

where T v is the contour in the +-sheet of Ei between and — p 2 along the real axis, and T w 
is the part of the real axis in the upper sheet of E 2 between — oo and — p 2 . The formula for 
u w shows that it does matter whether the equatorial plane is approached from the upper or 
the lower side (the Ernst potential is not regular at the disk). For p > p , we have u w = 



(G — 1 in the exterior of the disk). Similarly we get for the integral in the exponent of (|6.22| ) 



. \nGduj 00 + 00 - = ——7 f In Gd^oo+oo- . (7.29) 



27 



Summing up we can write the Ernst potential in the equatorial plane in the form 
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where $ = 1, and where v(P) = J P p 2 dv. The reality properties of the above theta functions 
imply together with ( |7.28|) the condition for equatorial symmetry /(— () = /(C)- Thus the 
imaginary part of / jumps at the disk. For p > p (where u w = 0), only the terms with even 
characteristics in (|7.30| ) will survive which leads to a real Ernst potential. This implies that 
the Ernst potential is regular outside the disk as it should be. The formula (|7.30|) can also 
be used to determine asymptotic quantities as angular momentum and ADM-mass in the 
limit of p — > oo as was done previously on the axis. 

A similar reduction as in the equatorial plane is possible on the axis. There the Riemann 
surface E' also has the involution T which makes it possible to factorize the surface into 
the surfaces E' x and E2 where the E< are as above and where E' x is Si with the cut [0, — p 2 ] 
removed. Thus the theta function 0' on the surface E' of genus 2s — 1 can be expressed via 
theta functions on surfaces of genus s — 1 and s respectively. In the case g = 2, this will not 
lower the genus of the Riemann surfaces under consideration on the axis. We will not give 
the Ernst potential on the axis since it will not be used here (the formula is helpful if one 
wants to calculate the multipole moments on the axis for higher genus). 



8 The case g = 2 



The simplest non-static solutions within the class (|6.22 ) are of genus 2 since solutions of genus 
belong to the Weyl-class. Interestingly the first solution to a physically relevant boundary 
value problem, the rigidly rotating dust disk |15[ with dust parameter v = 2Q 2 p^e~ 2Vo where 
Q is the angular velocity in the disk, and e~ v ° — 1 is the central redshift and radius p n , 



1 1 



belongs to this subclass. There the characteristic is 

E = yji/v - 1, and the function G has the form G = (^l + v 2 [K 2 + l) 2 + v{K 2 + 1)^ 
where we have used dimensionless coordinates p/p Q and C/Po- Therefore we will discuss the 
case g = 2 as an example in more detail. 

Solutions (|6.22j ) of genus 2 will be regular except at the contour T z if 0(u;(oo~) +u) ^ 0. 
This is equivalent to the condition that the divisor A, defined by the Jacobi inversion problem 
ou(A) —u(D) = u, does not contain both oo~ and X where X = Pq or X = Pq. This implies 
that the equations 
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1 r In Gdr 
27ri Jr p(r) 

1 r InGrdr 
27ri Jr p(r) 



(8.31) 
(8.32) 



must not hold simultaneously for all p and ( in the vacuum region. Since the limits of the 
integration are fixed, this inversion problem will in general not have a solution. The equations 
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constitute a relation between the physical parameters that characterize the solution. In the 
case of the dust disk, this is the parameter v and the radius pq. This implies that the above 
condition will determine the allowed parameter range for v for which there are no further 
singularities in the whole spacetime except the disk. 

The condition for ergospheres can be obtained from a similar system of equations as 
above: simply replace oo~ in ( |B.31| ) and ( |3.32j ) by a point on Lu which must not be oo~ or X 
but is otherwise arbitrary. Then an ergosphere occurs if both equations hold simultaneously 
which gives in the above example the values for v for which there exists a non-empty set 
of points p and (, the ergosphere. Since one of the points in the divisor A is in this case 
essentially arbitrary, the conditions for ergospheres will be satisfied much more frequently 
than the condition for a singularity where both points of A are prescribed. 

The ergosphere has only common points with the axis in the ultrarelativistic limit which 
is given by ^(i^) = (we use the notation for elliptic theta functions of p6| ). This condition 
is equivalent to 

1 r InGdr . pE dr , 

/ — r^ = ( 2n + 1 ) — rr 8-33 

2vri Jr p'(r) K J J-e p'(r) v ; 

where n = 0, 1, 2, Since this relation is independent of the physical coordinates, it deter- 
mines the values for v at which the ADM-mass diverges. Together with the conditions (|8.31|) 
and ( |S.32|) , this determines the allowed parameter range for v: the absence of singularities 



except the disk and the fact that the mass shall vary between and oo (the ultrarelativistic 
limit). 

In the case g = 2, the potential in the equatorial plane and on the axis can be expressed 
in terms of elliptic theta functions. The formula for the axis reads 

f(p = o,0 
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In the equatorial plane we have 



- _ ^ + uy] 2n 1 )^K; 2n 2 ) - $ 2 {v + Uv] zn^hi^ 2n 2 ) Iv 

1 1 ? 3 (««-'y;2n 1 )i?4(««,;2n 2 )+7? 2 (« 1> -v;2n 1 )i?i(« 11 ,;2n 2 ) e [ ' ' 

where v = J^°Ji dv. For p > p , the exterior of the disk, we get 

V 3 (u v - v; 2TI 1 ) 

In both cases the formulae for a different characteristic are obtained by complex conjugation. 

The above relations illustrate that important features of the solutions of genus 2 can be 
discussed with the help of the standard elliptic theory. It is thus an interesting question 
which boundary value problems lead to solutions within this subclass. 
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9 Outlook 



In this paper, it was shown that it is possible to identify within a class of hyperelliptic 
solutions to the Ernst equation a subclass whose solutions could describe the exterior of a 
body of revolution: they are asymptotically flat, equatorially symmetric and regular except 
at the surface of the body. This subclass could consequently be interesting in the context 
of boundary value problems for the Ernst equation. There the boundary data are either 
induced by an interior solution or a surfacelike matter distribution as in the case of the dust 
disk: in the general case, one would have to fix two real functions at the boundary in order 
to satisfy the boundary conditions. Within the class considered here, one has the freedom 
to choose one real valued function (G) and a set of free parameters E*, the branch points 
of the Riemann surface. Whether the subclass discussed here can actually be used to solve 
boundary value problems, and when these degrees of freedom will be indeed sufficient is 
an open question. However it is remarkable that it is possible to identify a whole generic 
subclass of regular equatorially symmetric solutions. This gives reasonable hope that further 
solutions to physically interesting boundary value problems may be found within the class 
of hyperelliptic solutions to the Ernst equation. 
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